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An interesting strategy for indirect detection of Dark Matter comes through the amounts of electrons and 
positrons usually emitted by DM pair annihilation. The e + e" gyrating in the galactic magnetic field then 
produce secondary synchrotron radiation. The radio emission from the galactic halo as well as from its expected 
substructures if compared with the measured diffuse radio background can provide constraints on the physics of 
WIMPs. In particular one gets the bound of (<tav) = 10~ 24 cm 3 s -1 for a DM mass m x = 100 GeV even though 
sensibly depending on the astrophysical uncertainties. 



1. Introduction 

Among the indirect DM detection channels, 
the radio emission due to secondary electrons or 
positrons can represent a chance to look for DM 
annihilation. During the process of thermaliza- 
tion in the galactic medium the high energy e + 
and e~ release secondary low energy radiation, in 
particular in the radio and X-ray band, which in 
principle could be detected. Furthermore, while 
the astrophysical uncertainties affecting this sig- 
nal are similar to the case of direct e + , e~ de- 
tection, the sensitivities are quite different, and, 
in particular, the radio band allows for a the dis- 
crimination of tiny signals even in a background 
many order of magnitudes more intense. 

Indirect detection of DM annihilation through 
secondary photons has received recently an in- 
creasing attention, exploring the expected sig- 
nature in X-rays [11213) . at radio wavelengths 
[4|5|6|7j , or both |8|9|10j . In this paper we will 
focus our analysis on the radio signal expected 
from the Milky Way (MW) halo and its substruc- 
tures. The results and the details of the following 
approach can be found in Ref. [IT] . 



2. Dark matter distribution 

Our knowledge of the DM spatial distribution 
on galactic and subgalactic scales has greatly im- 
proved thanks to recent high resolution zoomed 
N-body simulations [15116117118] . These sim- 
ulations indicate that for the radial profile of 
the galactic halo the usual Navarro-Frank- White 
(NFW) distribution [19] 
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still works as a good approximation over all the 
resolved scales. The NFW profile is in fair agree- 
ment with the observed Milky Way rotation curve 
[20 , although, depending on the employed model, 
it is possible to find an agreement for many dif- 
ferent DM profiles. We emphasize, however, that 
the various profiles differ mainly in the halo center 
(for r < 1 kpc) where the uncertainties, both in 
numerical simulations and from astrophysical ob- 
servations are maximal. Thus, our analysis which 
explicitly excludes the galactic center, does not 
crucially depend on the choice of the profile. 

The parameters describing the halo are then 
determined imposing the DM density to be equal 
to ps = 0.365 GeV c~ 2 cm~ 3 near the Solar Sys- 
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tem, at a galactoccntric distance of Rs = 8.5 kpc. 

Simulations, however, predict a DM distribu- 
tion sum of a smooth halo component, and of an 
additional clumpy one with total masses roughly 
of the same order of magnitude. For the subhalo 
population we will assume a mass distribution 
oc Ti~ ; 2 and that they arc spatially distributed 
following the NFW profile of the main halo. The 
mass spectrum number density of subhaloes, in 
galactoccntric coordinates r, is thus given by 
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where A is a dimensional normalization constant. 

By using the constraints described in one 
can fix the values of free parameters rn, Ph and A, 
hence obtaining ru = 14.0 kpc, ph = 0.572 GeV 

c -2 cm -3 and A = 1A Q x 1Q -19 kpc -3 M -l_ 

A further piece of information is required to 
derive the annihilation signal from the clumps, 
namely how the DM is distributed inside the 
clumps themselves. We will assume that each 
clump follows a NFW profile whose parameters 
are fixed assuming the clumps concentration ac- 
cording to [21] . 

3. Radio Data 

In the following we will derive constraints on 
the DM emission comparing the expected diffuse 
emission from the smooth halo and the unresolved 
population of clumps with all sky observation 
in the radio band. In the frequency range 100 
MHz-100 GHz where the DM synchrotron signal 
is expected, various astrophysical processes con- 
tribute to the observed diffuse emission. Com- 
peting synchrotron emission is given by Cosmic 
Ray electrons accelerated in supernovae shocks 
dominating the radio sky up to ~ 10 GHz. At 
higher frequencies the Cosmic Microwave Back- 
ground (CMB) and its anisotropies represent the 
main signal. However, thanks to the very sensi- 
tive multi-frequency survey by the WMAP satel- 
lite, this signal (which represents thus a back- 
ground for DM searches) can be modeled in a 
detailed way and can thus be removed from the 
observed radio galactic emission [52]. Other pro- 
cesses contributing in the 10-100 GHz range are 



given by thermal bremsstrahlung (free-free emis- 
sion) of electrons on the galactic ionized gas, and 
emission by small grains of vibrating or spinning 
dust. 

In the following our approach will be to com- 
pare the DM signal with the observed radio emis- 
sion where only the CMB is modeled and re- 
moved. For this purpose we use the code de- 
scribed in |23| where most of the radio survey ob- 
servations in the range 10 MHz-100 GHz are col- 
lected and a scheme to derive interpolated, CMB 
^2) cleaned sky maps at any frequency in this range 
is described. 

4. DM Synchrotron Signal 

In a standard scenario where WIMPs ex- 
perience a non exotic thermal history, a 
typical mass range for these particles is 
50 GeV < m x < 1 TcV, while a simple estimate 
for their (thermally averaged) annihilation cross 
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section yields (ctav) = 3 x 10 _27 cm 3 s~ 
giving (<jav) ~ 3x 10 _26 cm 3 s~ 1 for VL cdm h 2 k, 0.1 
as resulting from the latest WMAP measure- 
ments. However, this naive relation can fail badly 
if, for example, coannihilations play a role in the 
WIMP thermalization process [M], an d a much 
wider range of cross sections should be consid- 
ered viable. In this work we consider values of m x 
from about 10 GeV to about 1 TeV, and ((Jav) 
in the range 10~ 26 10 -21 cm 3 s _1 

The e + e~ annihilation spectrum, 
dN e /dE e (E e ), for a given super-symmetric 
WIMP candidate can be calculated for example 
with the DarkSUSY package [25]. Once injected 
in the galaxy, by neglecting diffusion (see Ref. 
[TT] for details), the emitted e + e _ follow the 
steady state distribution which reads 
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is the constant rate at which DM annihilation 
injects new electrons in the galaxy and r = 
E e /b(E e , r) is the cooling time, resulting from the 




sum of several energy loss processes that affect 
electrons. In the following we will consider syn- 
chrotron emission and Inverse Compton Scatter- 
ing (ICS) off the background photons (CMB and 
starlight) only, which are the faster processes and 
thus the ones really driving the electrons equilib- 
rium. Other processes, like synchrotron self ab- 
sorption, ICS off the synchrotron photons, e + e~ 
annihilation, Coulomb scattering over the galac- 
tic gas and bremsstrahlung are generally slower. 
They can become relevant for extremely intense 
magnetic field, possibly present in the inner par- 
sees of the galaxy [5], and thus will be neglected 
in this analysis. 

The synchrotron spectrum of an electron gy- 
rating in a magnetic field has its prominent peak 
at the resonance frequency 
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This implies that, in practice, a ^-approximation 
around the peaks works extremely well. Us- 
ing this frequency peak approximation, the syn- 
chrotron cmissivity can be defined as 

This quantity is then integrated along the line of 
sight to get the final synchrotron flux across the 
sky: 
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where (7, b) are coordinates on the sphere and s 
the line of sight coordinate. 

The two contributes mainly differ due to the 
fact that, while the halo emission is proportional 
to the square of the NFW profile, the clumps one 
is simply proportional to it. (see Ref. [TT]). 

The halo component dominates in the central 
region of the galaxy, inside a disk with a radius 
of about 30 degrees, while the clumps component 
represents the main contribution in the external 
region (see fig. EI). 

5. DM Annihilation constraints 

The pattern and intensity of the DM radio 
map resulting from the sum of the contributions 




Figure I. Sky map of the galactic radio signal 
generated by the DM smooth halo and unresolved 



clumps at the frequency of I GHz for 



100 



GeV and (a A v) = 3 x Kr 26 cm 3 s _1 . The pecu- 
liar shape of the signal, pinched approximately at 
±30° and ±60°, reflects basically the structure of 
the magnetic field as seen in projection from the 
Solar System, where the observer is located. 



from the smooth halo and unresolved clumps is 
shown in Fig. Q]for m x = 100 GeV and (uav) = 
3 x 10 _26 cm 3 s _1 . Similar maps are obtained at 
different frequencies and different m x and (uav) 
to obtain DM exclusion plots. For our analysis 
we use a small mask covering a I5°xI5° region 
around the galactic center where energy loss pro- 
cesses other than synchrotron and ICS start pos- 
sibly to be relevant. We include the galactic plane 
although this region has basically no influence for 
the constraints on the DM signal. 

In Fig [3] we show the radio constraints on the 
DM annihilation signal in the m x -(<7Av) plane for 
various frequencies and various choices of the fore- 
ground. Several comments are in order. First, we 
can see that, as expected, the use of the haze at 23 
GHz gives about one order of magnitude better 
constraints with respect to the synchrotron fore- 
grounds at the same frequency. However, using 
also the information at other frequencies almost 
the same constraints can be achieved. This infor- 
mation in particular is complementary giving bet- 
ter constraints at lower DM masses. This is eas- 
ily understood since a smaller DM mass increases 
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Figure 2. DM synchrotron profile for the halo and 
unresolved substructures and their sum at 1 GHz 
for m x = 100 GeV and (a A v) = 3 x 10~ 26 cnrV 1 . 
The astrophysical observed emission at the same 
frequency is also shown. The gray band indicates 
the angular region within which the DM signal 
from the host halo dominates over the signal from 
substructures modeled as in section [2j 



the annihilation signal (oc ra~ 2 ) at smaller ener- 
gies, and thus smaller synchrotron frequencies. In 
particular, the constraints improve of about one 
order of magnitude at m x ~ 100 GeV from 23 
GHz to 1 GHz while only a modest improvement 
is achieved considering further lower frequencies 
as 0.1 GHz. This saturation of the constraints 
is due to the frequency dependence of the DM 
signal, that below 1 GHz becomes flatter than 
the astrophysical backgrounds so that the frac- 
tion of contribution from DM is maximal at about 
1 GHz. 

For low masses the constraints come more and 
more from lower frequencies. For example for a 
WIMP of 30 GeV the data at 100 MHz are 2 
orders of magnitude more constraining than the 
data at 10 GHz. For a WIMP of 1 GeV, from 
Eqj5] with a magnetic field of 0(/iG) only fre- 
quencies < 10 MHz would be useful to place con- 
straints on the DM signal. Although observations 



at this frequency exist [23], in general the sur- 
vey sky coverage is quite incomplete and the data 
quality is non-optimal. Observations in this very 
low frequency range should substantially improve 
with the next generation radio arrays LOFAR and 
SKA. 

6. Summary and conclusions 

Using conservative assumptions for the DM dis- 
tribution in our galaxy we derive the expected 
secondary radiation due to synchrotron emission 
from high energy electrons produced in DM an- 
nihilation. The signal from single bright clumps 
offers only poor sensitivities because of diffusion 
effects which spread the electrons over large ar- 
eas diluting the radio signal. The diffuse signal 
from the halo and the unresolved clumps is in- 
stead relevant and can be compared to the radio 
astrophysical background to derive constraints on 
the DM mass and annihilation cross section. 

Constraints in the radio band, in particu- 
lar, are complementary to similar (less stringent 
but less model dependent) constraints in the X- 
ray/gamma band |26|27j and from neutrinos [28j . 
Radio data, in particular, are more sensitive in 
the GeV-TeV region while neutrinos provide more 
stringent bounds for very high DM masses ( > 10 
TeV). Gammas, instead, are more constraining 
for m x < 1 GeV. The combination of the various 
observations provides thus interesting constraints 
over a wide range of masses pushing the allowed 
window significantly near the thermal relic possi- 
bility. 

More into details, we obtain conservative con- 
straints at the level of (ctav) ~ 10~ 23 cm 3 s -1 for a 
DM mass m x = 100 GeV from the WMAP Haze 
at 23 GHz. However, depending on the astro- 
physical uncertainties, in particular on the as- 
sumption on the galactic magnetic field model, 
constraints as strong as (a A v) ~ 10 -25 cm 3 s _1 
can be achieved. Complementary to other works 
which employ the WMAP Haze at 23 GHz, we 
also use the information in a wide frequency band 
in the range 100 MHz- 100 GHz. Adding this 
information the constraints become of the or- 
der of {(tav) ~ 10 _24 cm 3 s _1 for a DM mass 
m x = 100 GeV. The multi-frequency approach 
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Figure 3. (Left) Constraints in the m x -{aAv) plane for various frequencies without assuming syn- 
chrotron foreground removal. (Right) Constraints from the WMAP 23 GHz foreground map and 23 
GHz foreground-cleaned residual map (the WMAP Haze) for the TT model of magnetic field (filled 
regions) and for a uniform 10 //G field (dashed lines). 



thus gives comparable constraints with respect 
to the WMAP Haze only, or generally better 
for m x < 100 GeV where the best sensitivity is 
achieved at ~ GHz frequencies. 

The derived constraints arc quite conservative 
because no attempt to model the astrophysical 
background is made differently from the case of 
the WMAP Haze. Indeed, the Haze residual map 
itself should be interpreted with some caution, 
given that the significance of the feature is at the 
moment still debated and complementary anal- 
yses from different groups (as the WMAP one) 
miss in finding a clear evidence of the feature. 
Definitely the multi-frequency approach will be 
necessary to test in a convincing way a possible 
DM signal like the claim related to the WMAP 
Haze. Progresses are expected with the forthcom- 
ing data at high frequencies from Planck and at 
low frequencies from LOFAR and, in a more dis- 
tant future, from SKA. These surveys will help 
in disentangling the various astrophysical contri- 
butions thus assessing the real significance of the 
Haze feature. Further, the low frequency data in 



particular, will help to improve our knowledge of 
the galactic magnetic field. Progresses in these 
fields will provide a major improvement for the 
interpretation of the DM-radio connection. 



Acknowledgments 

G.M. acknowledges supports by the Spanish 
MICINN (grants SAB2006-0171 and FPA2005- 
01269) and by INFN-I.S.Fa51 and PRIN 2006 
"Fisica Astroparticcllare: Neutrini ed Uni verso 
Primordiale" of Italian MIUR. 

REFERENCES 

1. L. Bergstrom, M. Fairbairn and L. Pieri, 
Phys. Rev. D 74 123515 (2006) 
|arXiv:astro-ph/0607327| . 

2. M. Regis and P. Ullio, Phys. Rev. D 78 (2008) 
043505 |arXiv:0802.0234l [hep-ph]]. 

3. T. E. Jeltema and S. Profumo, 
larXiv:0805.T054l [astro-ph]. 

4. P. Blasi, A. V. Olinto and C. Tyler, 



6 



Astropart. Phys. 1 8 (2003) 649 
|arXiv:astro-ph/0202049| . 

5. R. Aloisio, P. Bl asi and A. V. Olinto, JCA P 
0405 (2004) 007 [arXiv:astro-ph/0402588| . 

6. A. Tasitsiomi, J. M. Siegal-Gaskins and 
A. V. Olinto, Astropart. Phys. 21 (2004) 637 
|arXiv:astro-ph/0307375] . 

7. L. Zhang and G. Sigl. larXiv:0807.3429l [astro- 
ph]. 

8. S. Colafrancesco, S. Profumo and P. Ul- 
lio, Astron. Astrophys . 455 (2006) 21 
|arXiv:astro-ph/0507575] . 

9. S. Colafrancesco, S. Profumo and P. Ul- 
lio, Phys. Rev. D 75 (2007) 023513 
|arXiv:astro-ph/0607"073| . 

10. E. A. Baltz a nd L. Wai, Phys. Rev. D 70 
(2004) 023512 [arXiv:astro-ph/0403528| . 

11. E. Borriello, A. Cuoco and G. Miele, 
larXiv:0809.2~990l [astro-ph] . 

12. D. Hooper, D. P. Finkbeiner and 
G. Dobler, Phys . Rev. D 76 (2007) 083012 
|arXiv:0705.3655l [astro-ph]] . 

13. D. Hooper, Phys . Rev. D 77 (2008) 123523 
|arXiv:0801. 43781 [hep-ph]]. 

14. P. Grajek, G. Kane, D. J. Phalcn, A. Pierce 
and S. Watson, |arXiv:0807.1508l [hep-ph]]. 

15. J. Diemand, B. Moore and J. Stadel, Nature 
433 (2005) 389 |arXiv:astro-ph/0501589| . 

16. J. Diemand, M. Kuhlen and P. Madau, 
Astrophys. J. 65 7 (2007) 262 
|arXiv:astro-ph/0611370] . 

17. M. Kuhlen, J. Diemand and P. Madau, 
larXiv:0805.4416l [ astro-ph]. 

18. V. Springel et al. larXiv:0809.0898l [astro-ph] . 

19. J. F. Navarro, C. S. Frenk and 
S. D M. White, Astrophys . J. 490 (1997) 
493 |arXiv:astro-ph/9611107| . 

20. A. Klypin, H. Zhao and R. S. Somerville, 
Astrophys. J. 57 3 (2002) 597 
|arXiv: astro-ph/0 1 10390] . 

21. J. S. Bullock et al, Mon. Not. 
Roy. Astron. Soc. 321 (2001) 559 
|arXiv:astro- ph/9908159] . 

22. M. Tegmark, A. de Oliveira-Costa and 

A. Hamilton, Phys. Rev D 68 (2003) 123523 
|arXiv:astro-ph/0302496l . 

23. A. de Oliveira-Costa, M. Tegmark, 

B. M. Gaenslcr, J. Jonas, T. L. Landecker 



and P. Reich. larXiv:0802. 15251 [astro-ph] . 

24. K. Griest and D. Seckel, Phys. Rev. D 43 
(1991) 3191. 

25. P. Gondolo, J. Edsjo, P. Ullio, L. Bcrgstrom, 
M. Schelke and E. A. Baltz, JCA P 0407 
(2004) 008 [arXiv:astro-ph/0406204 . 

26. G. D. Mack, T. D. Jacques, J. F. Beacom, 
N. F. Bell and H. Yuksel, larXiv:0803.0~157l 
[astro-ph] . 

27. M. Kachelriess and P. D. Scrpico, Phys. Rev. 
D 76 (2007) 063516 [arXiv:0707.0209l [hep- 
ph]]. 

28. H. Yuksel, S. Horiuchi, J. F. Beacom and 
S. Ando, Phys. Rev. D 76 (2007) 123506 
|arXiv:0707.01~96l [astro-ph]] . 



